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ABSTRACT 

The quantitative analysis of low resolution spectra of A and B supergiants 
is used to determine a distance modulus of 24.99 ±0.10 mag (995 ±46 Kpc) 



^Based on VLT observations for ESO Large Programme 171.D-0004. 
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to the Local Group galaxy WLM. The analysis yields stellar effective tempera- 
tures and gravities, which provide a distance through the Flux weighted Gravity- 
Luminosity Relationship (fglr). Our distance is 0.07 mag larger than the most 
recent results based on Cepheids and the tip of the RGB. This difference is 
within the la overlap of the typical uncertainties quoted in these photomet- 
ric investigations. In addition, non-LTE spectral synthesis of the rich metal line 
spectra (mostly iron, chromium and titanium) of the A supergiants is carried out, 
which allows the determination of stellar metallicities. An average metallicity of 
-0.87 ±0.06 dex with respect to solar metallicity is found. 

Subject headings: galaxies: — distances galaxies: abundances — galaxies: indi- 
vidual (WLM) — stars: early-type, supergiants 



1. Introduction 



WLM is one of the faintest dwarf irregular galaxies {Mb = —14), located in an isolated 
part of the Local Group. Detailed photometric studies have shown that it consists of a 



young population concentrated in a disk and an old extended metal poor halo (jFerraro et al. 



19891 : iMinniti fc Ziilstralll997l : lDolphinll2000l : iReikuba et al.ll200d : iMcConnachie et al.ll2005r ). 

The color of the red giant branch of the old population indicates a metallicity of [Fe/H] = 
-1.45 dejfl representing the end of the first star formation episode. On the other hand, the 
metallicity of the you ng population is soniewhat ambiguous. While measuremen ts from a 
number of Hll regions (jSkillman et al.lll989l : iHodge &: Millerlll995l : iLee et al.ll2005l ) yielded a 
low nebular oxygen abundance o f [0/H] =-0 . 8 dex , a detailed high resolution spectroscopic 
study of an A-type supergiant (jVenn et al. 2003 ) obtained [0/H] =-0.2 dex indicating a 
large discrepancy be tween stellar and nebular chemical composition. On the other hand, 
Bresolin et al.l (120061 ) in their low resolution (5 A) spectroscopic survey for supergiants in 
WLM studied three early B supergiants and found an average value of [0/H] = -0.8 dex 
comparable with the Hll regions. 



The spectroscopic sample gathered in the survey by iBresolin et al.l (120061 ) also con- 
tained high signal-to-noise data of several late B and A supergiants (B5-A7) of luminosity 
class la and II. These objects with their rich spectra of metal lines (mostly iron, titanium, 
chromium) are ideal for a further investigation of stellar metallicity. However, at the time 



^We use the common notation [X/Y] = log(X/Y) - log(X/Y)Q, with the solar abundanc e s giv en by 
Grevesse fc Sauvall (|l998l ). except for oxygen, for which we adopt the value from lAsplund et al.1 (|2004l ). 
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of the publication of this work the density of metal lines in the optical spectra together 
with the low spectral resolution did not allow for a determination of stellar parameters using 
the standard analysis techniques appli e d for the early-type B supergiants. Recently, this 
situation has changed. iKudritzki et al.l (l2008l . hereafter K08) in their study of the Sculptor 
spiral galaxy NGC300 at a distance of 1.9 Mpc developed a new technique to quantita- 
tively analyze low resolution spectra of A supergiants with respect to stellar parameters and 
metallicity. This technique makes use of the Balmer jump and the Balmer lines to obtain 
effective temperature and gravity and determines metallicity from a x^-minimization of ob- 
served and theo retical spectra i n selec ted spectral windows. It can now be applied to the data 
set obtained by iBresolin et al.l (120061 ) and allows for a more comprehensive straightforward 
spectroscopic study of the metallicity of the young stellar population in WLM. 

The distance to WLM has been determined from photometric studies of the old stellar 
population using the tip of the red giant branch (TRGB), the horizontal branch (HB) and 
full color-rnagnit ude diagrams, yielding distance moduli between 24.7 and 24.95 mag (see 
Rizzi et al.l 12007, and referen c es ab ove). The most recent comprehensive multi-color survey 
for Cepheids by lGieren et al.l (l2008l ). including J- and K-band photometry, found a distance 
modulus of 24.92 mag with a very small random error of 0.04 mag and a systematic error of 
0.05 mag. These authors determined an average interstellar reddening E(B-V)=0.08 mag, 
significantly larger than the Galactic foreground reddening and, thus, relevant for the derived 
distance. Under such circumstances, an independent distance and reddening determination 
using the young stellar population is ideal for verifying photometrically based distances. 

Such an independent distance determination c an be provided t hroug h the quantita- 
tive spectral analysis of the B and A supergiants. IKudritzki et al.l (120031 ) and K08 have 
shown that a tight correlation exists between the absolute bolometric magnitude and the 
flux weighted gravity g/T^^f. This relationship is predicted by stellar evolution theory and 
allows distance determinations with a precision comparable to Cepheids. 

In this paper, we will carry out a spectral an alysis of the low resol ution optical spectra 
of 8 late B and A supergiants (B5-A7) observed by lBresolin et al.l (120061 ) to determine stellar 
temperatures, gravities, luminosities, masses and metallicities. For the analysis we will use 
the basic concept introduced by K08, however, in a modified form. We will present a new 
numerical fit algorithm based on an empirical application of the Karhunen-Loeve expansion 
usually referred to as Principal Component Analysis, which allows for an automated analysis 
of a large number of objects using large comprehensive grids of model atmos phere spectra. 
We will then combine the results with those obtained by lBresolin et al.l (120061 ) for the early- 
type B supergiants and determine an independent distance using the flux-weighted gravity 
luminosity relationship, fglr. This will be the first distance determination using this new 
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distance determination method. 

The paper is organized as follows. After a brief description of the observations in §5] we 
present the modified analysis method in §31 In §11 we test our method by comparing with 
the results obtained in a previous detailed high resolution study for one of our targets. The 
results of the analysis, i.e. the derived Tejj, log g, metallicity, E(B-V), radius, luminosity and 
mass, are given in §51 The distance determination using the fglr is carried out in §6] and 
§7] will present a discussion of all the results. 



2. Observations 



The spectroscopic observ ations were obtained wit h the Focal Reducer and Low Disper- 
sion Spectrograph 2 (F0RS2. lAppenzeller et al.lll998l ) at the ESO Very Large Telescope in 
multi-object spectroscopy mode in one night of good seeing conditions (better than 0.7 arc- 
sec) on 2003 July 28 as part of the Auracaria Project (iGieren et al.ll2005l ). The total exposure 
time is 4500s and the air mass during the observations was smaller than 1.07. The spectral 
resolution is approximately 5 A and the spectral coverage extends over 2500 A centered 
at 4500 A for most of our targets. The spectra were also flux calibrated so that spectral 
energy distributions, in particular the Balmer jump, can be u sed to constrai r i the stellar 
parameters. The observational data set has been described by iBresolin et al.l (120061 ). The 
paper also contains finding charts, coordin ates, photometry (V, I from the Las Campanas 
1.3m telescope, see iPietrzyhski et al.l 120071 ) . radial velocities and spectral types. Of the 19 
confirmed supergiants found with spectral types ranging from late O to G, 8 were of spectral 
type B5 to A7. The S/N per pixel for these objects is between 40 to 120, sufficient for our 
quantitative analysis. This sample listed in Table [J an d selec ted for this work. The target 
designation follows the nomenclature in iBresolin et al.l (120061 ). Note that we consider only 
stars of their A set because of the significantly better S/N. 

Along with the broad band photom etry in the afo r emeii tioned reference, we also use 
B-, V-, R- and I-band data published by iMassey et al.l (120071 ) to constrairi the in terstellar 
reddening E(B-V), as well as J- and K-band photometry from lGieren et al.l (120081 ). The IR 
photometry, available only for some stars, is compiled in Table [2l 

Prior to the quantitative spectral analysis, we carefully explore the dataset available for 
each star, to identify spectral regions showing signs of contamination due to imperfect sky 
subtraction, nearby cosmic rays, and similar other observational effects. These regions are 
masked out and not used in the analysis. 
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3. Quantitative spectral analysis 

As explained in K08 the principal difficulty in the analysis of low resolution spectra of A 
supergiants lies in the simultaneous determination of effective temperature and metallicity. 
The classic method based on the use of ionization equilibria does not work, since the required 
lines from the neutral species (Mg I, Fe I, . . . ) are in general very weak and disappear at low 
resolution within noise and the blends of spectral lines. While one could use the information 
from the stronger lines, which define the spectral type, and a calibration of spectral type 
with effective temperature, such a relationship is metallicity dependent, which requires the 
simultaneous and independent determination of metallicity. 

K08 were able to solve this problem by using the Balmer jump as an independent 
temperature indicator in addition to the Balmer lines, which can be used as a measure of 
gravity, and to the rich metal line spectrum, which yielded the metallicity through a x^- 
minimization comparing selected windows of the observed spectra with very sophisticated 
non-LTE line formation calculations. 

The methodology employed in the present work follows closely the ideas presented by 
K08, however with some significant modifications. While K08 constructed fit diagrams of 
the Balmer jump and the Balmer lines in the T^ff-log g plane through "by eye" fits of the 
observed line proffies and the spectral energy distribution (SED) and then determined the 
metallicity for the effective temperature and gravity at the intersection of the two fit curves. 
Here we have developed a completely automated numerical method, which finds the best fit 
of the model atmosphere synthetic spectra to the observed spectra. The advantages of this 
new method are obvious. It provides a more quantitative and objective way to obtain the 
fit and the automated procedure allows a quantitative analysis of a large number of spectra, 
which in the era of efficient multi-object spectrograph facilities is a very important aspect. 
We will describe the method in the subsections below. 



3.1. Brief description of the model grid 



The basis for the quantitative spectral analysis is the same comprehensive grid of line 
blanketed LTE model atmospheres and very detailed non-LTE line formation calculations 
used by K08. For a detai led description of the physics of the model atmosphere and line 
formation calculations see iPrzybilla et al.l (120061 ) . 



Our grid of models covers a range of 8300 < T^ff < 15000 K (with steps of 250 K and 
500 K below and above 10^ K, respectively) and 0.75 < log g < 2.60 (with steps of 0.05 dex) 
in effective temperature and gravity with the following metallicities at each grid point, [Z] 
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= log(Z/Z)Q: -1.30, -1.00, -0.85, -0.70, -0.60, -0.50, -0.40, -0.30, -0.15, 0.00, 0.15, 0.30. The 
quantity Z/Zf^ is the met allicity relative to the Sun. The s olar abundances were taken from 
Grevesse fc Sauvall (Il998l ). except for oxygen, which is from lAsplund et al.l (12004 ) . Based on 
trend in the literature, a microturbulence of 8 km s~^ is adopted for the low gravity models 
gradually changing to 4 km s~^ at higher gravity. For further details we refer to K08. 

This grid was primarily designed for the analysis of la and lb luminosit y class objects. 



Severa l of the WLM stars selected for the present work were classified by iBresolin et al. 



(I2OO6I ) as luminosity class II. To consider these objects, we enlarged the grid with two new 
effective temperatures, 8100 and 7900 K, and extended to higher gravities all the models 
with Teff < 10^ K [log g < 2.75 dex). 



3.2. General methodology 

The goal of the analysis is to determine for each star a set of four parameters {T^jf, logg, 
[Z] and E(B-V)}, for which the corresponding synthetic model spectrum best matches the 
observations. We have three data sets available: 1) normalized and rectified spectra, which 
will give line profile information for the Balmer and metal lines; 2) flux calibrated spectra 
providing information about the Balmer jump, and 3) broad band photometry which can be 
used for the construction of the observed long wavelength SED. 

K08 have shown in detail how the theoretical spectral line profiles, the SED and the 
Balmer jump depend on temperature, gravity and metallicity. Generally speaking, the 
Balmer lines depend mostly on gravity, but show also a significant temperature depen- 
dence, whereas the Balmer jump depends mostly on temperature, but also on gravity. The 
dependence on metallicity of the Balmer lines and the Balmer jump is very weak. Using 
these properties, two different curves can be defined in the T^ff-log g plane, one being the 
locus of the models that reproduce the Balmer lines when the temperature is adopted, and, 
conversely, another tracing the locus of the models that reproduce the Balmer jump once 
the gravity is assumed. The intersection of these two curves yields the gravity and the 
temperature of the star observed. The classic method used by K08 is to construct these 
curves through a by-eye fit of the models to the observed data. In our new approach we have 
decided to use an automated numerical fit of the model spectra to the observed data. 

a) Our algorithm is iterative and begins with the fit of the Balmer lines. We start with 
a first value of Tgff and [Z] and try to find the model gravity log g which best matches the 
shape of the observed Balmer profiles. For this purpose, we carry out a Principal Component 
Analysis (PGA) of the normalized model spectra in spectral windows around the Balmer 
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lines. For a fixed Te// and [Z] we have n models with different gravities at each grid point (n 
varies from 20 to 32 depending on the grid point), each of them providing a set of normahzed 
fiux values Sj(A) (where the index j runs over the n models) around a selected Balmer line. 
These model fiux values define a matrix A (with elements Ajj =Sj(A), where i labels the 
wavelength points and j runs over the models with different log g). Following the concept 
of PCA, we det e rmine the eigen values and the eigenvectors of the covariance matrix cov(A) 



( iDeemingi Il964j : IWhitneyi 119831 ) and identify the largest eigenvalue and its corresponding 
eigenvector as the principal component associated with the parameter log g. Fig. [T] shows 
the typical eigenvectors v(A) for three Balmer lines. As a consistency test we check that the 
cumulative percentage variance associated with this principal component is always above 
95% in all cases, meaning that more than 95% of the information in the input data can be 
described by the principal component alone. 

The projection of the model fiux matrix A onto the eigenvector v allows us to define a 
relationship between the log g values of the models and their projections ^{log g), for each 
Balmer line 

S, (A) V (A) = $, = $ {log g) 

A 

which are also shown in Fig. [H Projecting the observed normalized fiux values F\ in each 
Balmer line spectral window onto the corresponding eigenvector yields an observed value, 
which can then be compared to ^{log g) to find the best gravity log g at the selected values 
of Teff and [Z]. Comparing the values obtained from different Balmer lines (we usually use 
Hio to H/j) allows us to assign a mean value and an uncertainty u to this determination of 
the best fitting gravity. 

b) Adopting this pair of iTeff) log g) we can now determine a new metallicity by defining 
spectral windows with many metal lines. Fig. [2] shows a typical example of different spectral 
windows and their corresponding PCA eigenvectors, this time determined with respect to 
metallicity [Z]. The projections $([^]) and the observed value are also shown. Using different 
spectral windows we can again assign an average value of [Z] and an uncertainty lv. For this 
new [Z] we iterate the determination of log g, but because of the weak dependence of the 
Balmer lines on metallicity changes are usually small. 

c) In the previous steps we used Balmer and metal line profiles to obtain log g and [Z] 
(and their corresponding uncertainties) at a selected value of T^ff. Now we adopt these new 
values of log g and [Z] and use the observed SED around the Balmer jump to obtain a new 
value of Teff. This is again done by PCA (but see below) as shown in Fig. [31 The model 
atmosphere SED corresponding to this Tg//, log g and [Z] is used to determine the reddening 
E(B-V) from a comparison ( y ^-minimization ) with the observed broad band photometry. 



adopting the extinction law by ICardelli et al.l (Il989l ). as well as Ry = Ay / E{B ~V) = 3.1 
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Since the region used to determine Tg// could be affected by this extinction, we alternatively 
iterate the values of T^ff and E(B-V) until convergence in both is reached, defining the best 
Teff and the best E(B-V). Unlike the cases of surface gravity and metaUicity, the temperature 
is determined from a single feature, the Balmer jump. The uncertainty in Te/j is thus given 
by our ability to effectively measure the Balmer jump, and not from temperature values 
derived from different spectral windows. In order to estimate this uncertainty, we proceed 
in the following way. Two values of the temperature are computed simultaneously, by using 
two different definitions of the Balmer jump (see Fig. [3]). The first one (Fig. corresponds 
to the value derived from the application of the PCA base vector (shown in Fig. [3t), while 
the second value (Fig. [8)3) is based on the Balmer jump index Db introduced by K08, 



In the case of ""^ , N is the number of wavelength points considered for the mean. Note 
that in both cases, we are considering all the models in the grid for a given pair {log g, [Z]), 
as previously explained. In an ideal case, with infinite S/N, both values would be exactly 
the same. Once noise is present, there is a difference between both temperatures, that 
reflects how precisely the Balmer jump can be measured, defining therefore the uncertainty 
in temperature, ATg//. 

Apart from this ATg//, the difference between the initially assumed temperature and 
the final derived value defines an error 6, which can be assigned to the set of parameters 
{Teff,log g, [Z],E {B — V)}. This S can be interpreted as the distance between the guess 
temperature and the true T^jf, because of the weak dependence of the Balmer jump on 
gravity and metaUicity. 

We repeat the whole procedure (a, b, c) for different initial values of Te// and, thus, ob- 
tain a set of parameters vTj = {T^ff, log g, [Z], E (B — V)}j together with their corresponding 
weights {6i,uji). The final solution 11 = {T^ff, logg, [Z] and E(B-V)} is then obtained as a 
weighted mean 



where the indexes i,j run over all the individual solutions for the corresponding parameter 
{Teff, logg, [Z] or E(B-V)). In a similar way, a formal error (uncertainty) is obtained from 



Db 



(log(Fr^)) - (log(Ff-*)), where 

(log(Ff = {Eii \og{F^J}/N, 3585A < A, < 3627A , and 




n 



(A nf 



n)2 
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This analysis method can only be applied in those cases for which the Balmer jump 
region is observed. In some instances, due to the location of the star in the observed field and 
its spectrum on the detector, this area is not observed, and we have to rely on the normalized 
spectrum to obtain information about the effective temperature. In this situation as already 
discussed by K08, the relative strength of lines from different species can be used. If the star 
is hot enough to show He I lines (earlier than around AO), these can be used to constrain the 
temperature. The accuracy is usually comparable to that when the Balmer jump is used. 
If the star is cool enough (later than B9), the relative behavior of different metal lines (for 
example Ti II versus Fe II lines, or the strength of the Fe-Mg blend located at ~5150 A) 
could be used. To find the temperature, we then proceed in a similar way, defining several 
spectral windows in which these metal or helium lines are present, to construct the PCA 
eigenvectors (from the normalized synthetic spectra), but now with respect to temperature. 
There is an underlying assumption, though, when using this method: the adopted metallicity 
pattern. The implications of our adopted solar abundance pattern are discussed in Sect 13.41 
Temperatures determined by this method have a somewhat larger uncertainty, related to the 
accuracy with which the metallicity can be determined. 

Metallicity and microturbulence are closely related with respect to the values derived 
from any analysis. As discussed by K08, it is difficult to constraint microturbulence at such 
low spectral resolution, since lines from a given species, with very different strengths, are 
required. In consequence, as explained in §3.11 for each object, the microturbulence is given 
by the final {Teff,log g) pair. We will discuss in §3.41 the possible effects of the adopted 
relationship between the microturbulence and the {Teff,log g) values. 

3.2.1. Uncertainties 

The formal accuracy of the final solution depends primarily on the signal-to-noise ratio. 
Our experience shows that, for the low spectral resolution considered in this work, a minimum 
S/N~50 is required for the metallicity determination in order to keep the uncertainties in 
a tolerable range (see next section), while Teg and log g determinations, both based on 
stronger features (the Balmer jump and the Balmer lines) can still be carried out at much 
lower, S/N~15. 

The formal errors arising from the application of our algorithm, as previously defined, 
are rather small. Typically, surface gravity and metallicity can be constrained to better than 
0.05 dex, while the effective temperature uncertainty results in a few K. To define more 
meaningful errors, we also consider the uncertainties derived for the model in the grid with 
parameters closest to the final solution, i.e., we use the oj values (defined above) obtained 
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for the set {Te/fjog g, [Z]} closest to the final solution. These provide an idea about how 
well the observations are reproduced by a single model. 

The final errors quoted are obtained as a quadratic combination of these two sets of 
uncertainties. Typical values for these final errors result in ~2-5% in temperature, ~0.1- 
0.2 dex in surface gravity and ~ 0.2-0.3 dex in metallicity. As an example, Fig. H] shows the 
solution obtained for one of the stars in the sample, A14. The small box encloses the solution 
with the formal errors in the Te//-logg plane. As can be seen this solution area agrees with 
the intersection of the conventional fit curves for the Balmer jump and the Balmer lines, 
with the difference of the solutions well below the resolution limit of both methods. The 
final adopted errors are represented by the dashed box. 

With regard to the uncertainties in the fiux weighted gravity, it must be noted that the 
errors in log g and T^ff are correlated (see Fig. H]), which reduces the uncertainties in log gp 
in most cases. The reader is referred to K08 for a detailed discussion on this topic. 



3.2.2. Goodness- of- fit assessment 



With the final set of parameters known for each star, it is useful to evaluate the fitness 
of each final solution, with the goal of identifying problems not detected in the analysis. 
Given a pair observation-final model {0\, Mx), we define the residuals of a spectral window 
with rii wavelength points as 



r^{Mx,Ox) 



j=i 



1 1/2 



The fitness R (11) of a parameter 11 is given by the sum of these relative residuals over 
all the spectral windows considered in the determination of the corresponding parameter (as 
described above), weighted by the (normalized) S/N of each individual window 



i?(n) = ^r^Si 



, with 



In the case of the surface gravity and metallicity, the normalized spectra are used to calculate 
the fitness. For the effective temperature, the flux calibrated data and the synthetic SEDs 
are considered to evaluate the fitness when the Balmer jump is observed. If it is not observed 
the normalized spectra are used. 
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Alternatively, we could define the fitness for each window as 

where e is the error allowed in the fitting (tolerance), defined as e = S/N~^. The global 
fitness, for each parameter individually, is given by the sum over all the windows, in this 
case without weighting by the S/N since its effect is already accounted for with the tolerance. 
This second fitness evaluator has the advantage of containing additional information in its 
sign, with a positive/negative value reflecting whether the corresponding parameter is over 
or under-estimated. This fitness definition is particularly useful to detect variations in the 
goodness-of-fit from star to star. 

We define the global fitness of the model with respect to the observation, R (Ma, Ox), as 
the sum of the individual fitness values of each parameter, R (11). From the definitions, it is 
clear that the smaller the absolute value, the fitter the model. Fitness values corresponding 
to the final solutions are presented in Table [3l and they will be used in the following sections 
to discuss the results. 



3.3. Consistency tests 

We performed a number of simple tests to check the analysis algorithm. In all the 
following cases, we emulated the observed spectra by degrading the models to a resolution 
of 5 A (FWHM) and re-sampling them to 1.32 A (the spectral resolution and dispersion 
provided by F0RS2 when equipped with the 600B grism). Note that these tests are meant 
to check our ability to reproduce known input parameters. 

First, we verified that, for any given model in the grid, we were able to recover almost 
exactly its parameters. In a second step, we created a number of models with parameters 
within the limits of the grid, but without any corresponding model in the grid. Three 
different sets of Tg/j-logg, representatives of a late B, an early A and a mid A supergiants 
were considered, for four different metallicities. In all cases, we were able to recover the 
parameters to better than 1% in effective temperature (based on the Balmer jump), 0.02 dex 
in surface gravity and 0.03 dex in metallicity. These two tests were carried out without 
considering noise, since they were intended to probe the ability of our algorithm to recover 
known parameters, without any influence from external sources. 

A third set of tests was performed to check the effect of noise on the determination of 
metallicity, and hence to evaluate the minimum S/N required for (meaningful) metallicity 
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determinations. We selected a number of models in the grid representative again of the 
different spectral groups, degraded them to different S/N (100, 50, 30 and 15) and carried 
out a metallicity analysis. For each model and each S/N value, we did 100 independent trials. 
The results for the early A-type case are presented in Fig. O (the results are very similar for 
all the other Te//-logg pairs considered). The first two rows correspond to [Z] = 0.00 dex, 
while the other two rows are for [Z] = —0.85 dex. Each single plot displays the relative 
frequency of occurrence versus the difference between the input and recovered metallicities, 
presenting also the sigma of the distribution. Since this sigma only describes the dispersion 
introduced by the noise, the global uncertainty obtained in a generic analysis would be larger, 
once the effects of the uncertainties in the other parameters (temperature and gravity) are 
accounted for. As a general recipe, a maximum as/N ~ 0.10 dex at a fixed pair of Tg// and 
logg would be required in order to have a final global uncertainty cr ~ 0.20 dex. 

Clearly, the minimum S/N depends on the metallicity: the higher the [Z], the stronger 
the lines, hence the lower the S/N required to detect them. At the same time, it also depends 
to some extent on the spectral type; for a given metallicity, cool mid A-types (Al to A3) 
present stronger metal features (as well as a higher line density) than late B-/early A-types 
(B5 to AO). From Fig. [S], it is possible to identify a minimum S/N~50, in particular at low 
metallicities, to meet the requirement of a global metallicity uncertainty around 0.2 dex. 

3.4. Dependence of the results on hidden parameters 

In order to keep the task of creating such a huge model grid manageable, only three 
parameters (the most important ones, Tg//, logg and [Z]) are explored so far. In the following 
we want to discuss some of the possible implications regarding how the models are calculated. 

As previously stated, metallicity and microturbulence are coupled in any analysis. 
Aimed at exploring the effect that our assumption about microturbulence ^ (given for a 
pair Te//-logg) could have on the results, we calculated a small set of models with different 
values of ^ (but the same Tg//, logg, [Z] considered for the grid), and analyzed them in the 
same manner as the observations. For a variation of ^ of ±2 km s~^ (which is the typical 
dispersion obtained in high resolution analyses of these type of objects, for a given luminosity 
class) we found an error of less than 1%, 0.01 dex and 0.07 dex for effective temperature, 
surface gravity and metallicity, respectively. As expected, both Tgjj and logg are unaffected 
by a change in microturbulence, and the larger effect is produced in the metallicity. Here 
we want to stress again that we are not dealing with individual lines, which certainly could 
present large variations with S, (like for example Mg II 4481 A), but with all the spectral 
lines at once (or with as many of them as possible); the use of the whole spectrum weights 



-13- 



the effect of the turbulence differently, reducing its impact on the derived global metallicity. 
This is, on the other hand, a warning sign for the use of small sections of the spectrum 
(at low resolution) to determine individual elemental abundances, which would be seriously 
hampered by the unknown ^. 

We proceeded in a similar way in order to evaluate the impact of our assumption of a 
solar abundance pattern. We calculated a set of models for which the ratio of the ^-elements 
to Fe-group elements relative to the solar reference vary, with [a/Fe] = ±0.15, ±0.30 dex 
(note that by definition, the value used in the grid is [a/Fe] = 0.0 dex). The maximum change 
in effective temperature (based on the Balmcr jump) and surface gravity due to these changes 
in relative abundances is less than 1% and 0.01 dex respectively. The difference goes up to 
6% in temperature when using the normalized spectrum instead of the Balmer jump. This 
result is not surprising since in this case, T^ff is based on the relative strength of different 
species, in particular of Ti II to Fe 11 lines, which depends on the relative abundances as well 
as on the temperature. As a byproduct, this exploratory study allowed us to refine the limits 
for some of the spectral windows considered in the determination of metallicities, allowing 
also to identify features that arc produced mainly by Fe-group elements (Fe, Ni, Cr, . . . ) 
and, conversely, features produced by a-elements (in particular Ti). This will be used in the 
future to evaluate the possibility of extracting information from the spectra not only for the 
global metallicity but also for the relative abundances of some prominent species. 



3.5. Comments on individual stars 

In this section we present a brief discussion of the analysis of some of the targets. 

• A4: the Balmer jump region was not observed for this star, therefore we have to use 
the normalized spectrum to estimate T^ff from metal lines (the star is too cool to 
present He lines). In this situation, T^ff is coupled to [Z], the well know problem of the 
spectral type-metallicity degeneracy. Due to the combination of the low metallicity 
and the relatively low S/N, the errors in [Z] and T^ff are relatively large. 

• A5 and A17: the Balmer jump is missing for both stars. Fortunately, they are hot 
enough to present He I lines, and can be used to constraint Tg/j. Note that both 
stars have S/N above 50, which greatly improves the situation with respect to the two 
previous cases. 

• A6: the analysis of this star is problematic. While the determination of T^g and log g 
using the Balmer jump and the Balmer lines is straightforward and results in fitness 



-14- 



values similar as for the other stars, the determination of [Z] is not. The metalhcity is 
not well defined, as is indicated by the high values of ([Z]) and Q ([Z]) compared to 
the other stars. Moreo ver, the high va l ue of [ Z]=-0.15 dex is puzzling. We note that the 



spectral type found by iBresolin et al.l (120061 ) based on the relative strength of Ca II H 



and K lines is A7, suggesting a significantly cooler temperature than the one obtained 
from the Balmer jump. While there are other examples of A supergiants where the Ca II 
related spectral types do not match Ti,ff (see discussion of star A16) mostly because of 
the contribution by interstellar lines (note that E(B-V) of A6 is 0.29 mag), we have 
tried to obtain a spectral fit at cooler temperatures (Tg// = 8100 K, log g = 1.50 dex) 
with a metallicity [Z]=-0.7 dex more in line with the other stars. At these parameters, 
the Balmer lines are reproduced well, but the fit of the Balmer jump is unacceptable 
as indicated by the fitness values in Table [31 In addition, the fitnesses i?([Z]) and 
Q ([^]) are only marginally improved and still significantly worse than the fits for the 
other stars. The value for E(B-V) at this cool temperature is 0.27 mag, only slightly 
lower than before. 

The problem described above might be caused by stellar multiplicity. We have tried to 
emulate the observables (Balmer jump, normalized spectrum and broad band photom- 
etry) by combining two models, one hot to reproduce the Balmer jump and one cool to 
produce the rich metal spectrum and the Balmer lines. However, it has been impossible 
to find a combination of models to fit everything. Also, the star does not show signs 
of companions in archival WFPC2/HST images. Even more puzzling is the fact that 
all the different photometric measurements considered, from B- to K- band (including 



WFP C2/HST photometry in filters F555W and F814W taken from iHoltzman et al. 



20061 ) are satisfactorily reproduced with a single model solution, invoking high redden- 
ing. 

Another possibility is that we are dealing with a completely different type of star with 
just the photometric and (partial) spectroscopic appearance of an A supergiant. The 
only solution that comes to our min d is a low mass st ar, for instance in post-RGB 



or post-AGB phase, in the Galaxy. I Venn et al.l (120031 ) have already discussed this 



possibility and found that this is rather unlikely. Indeed, if we assume the typical 
mass of 0.5 M© for these objects, then with log g = 1.95 dex and T^jf = 8750 K, 
we obtain log L/Lq = 2.9 or Mboi = —2.5 mag, with which, an apparent magnitude 
of mboi = 18.9 mag would put the object at a distance of 190 Kpc, certainly outside 
the Galaxy. We also note that A6 has a radial velocity typical for a WLM member 



(IBresolin et all 120061 ). 

At this point, we consider the metallicity of A6 obtained in our analysis as unreliable. 
We will include the object in our discussion of the fglr and of interstellar reddening. 
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but will regard it uncertain. 



A14: this star has been previously studi ed by I Venn et al.l (20031 ). to which these authors 
referred as WLM-15. In their analysis, I Venn et al.l (120031 ) found a significant discrep- 
ancy between the oxygen abundance of the star and the abundance s obtained from 
H II regions. We compare our results with those of I Venn et al.l (120031 ) in the following 
section. Unfortunately, our wavelength coverage does not include O lines, therefore it 
is not possible to carry out an independent analysis of the oxygen abundance. 



A16: the spectral type classification by lBresolin et al.l (120061 ) . based on the Ca II lines, 
is not consistent with the spectrum, nor with the results of the analysis (see below). 
There is at least another well known case for which this also happens. The Galactic 
star HD 12953, the standard Al la in the MK system, would be classified as an A3 
based on the strength of its Ca II lines, as already discussed by lEvans fc Howarth 
( 120031 ). Unlike the case of A6, a completely consistent solution can be found. 



4. Comparison with high resolution spectroscopy 



One of the stars in our sample, A14, was previously studied by lVenn et al.l (120031 . here- 
after V03). These authors obtained a high resolution optical spectrum with UVES at the 
VLT (R~32000), and subsequently analyzed it using model atmosphere techniques similar 
to the ones employed here. This offers the possibility of a comparison with our low resolution 
work. 

For the fundamental parameters, T^ff and logg, V03 obtained 8300 ± 200 K and 1.60 ± 0.10 
dex, while our solution is 8270 ±145 K and 1.60 ±0.12 dex. The agreement is very good. 
Worth noticing is the fact that the effective temperature of V03 is based on the Mg ionization 
equilibrium, while ours is based on the Balmer jump. With respect to the chemical abun- 
dances, the high resolution work is far superior in the sense that it can provide individual 
elemental abundances, while we derive a global metallicity. In order to compare with our 
result, we used V03's individual abundances of Fe, Ti, Sc and Cr (and our solar references) 
and calculated a mean [Z] value of -0.44 ±0.13 dex (the latter number is the standard devi- 
ation of this weighted mean). Note that neither O nor N lines available to V03 are covered 
in our observed spectral range, thus we decided not to include them to compute [Z] . Within 
the uncertainties, this value is in good agreement with our derived value of -0.50 ± 0.19 dex. 
In fact, this good concordance can be seen in Fig. [HI were we show the UVES spectrum of 
A14 and a model with the parameters resulting from our low resolution analysis, convolved 
with the appropriate UVES instrumental profile and rotational velocity. 
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While the comparison can only be done for one star, the close agreement found is 
certainly very encouraging, in particular given the relatively low metallicity of the obj ect. 
This agreement is also in consonance with the result presented by iBresolin et al.l (120061 ) for 
the early B-type supergiant A9 in the same galaxy. 



5. Results 

The stellar parameters obtained in our work are summarized in Table |H In the following 
we discuss metallicity, interstellar extinction and stellar parameters. 



5.1. Metallicities 



The metallicity of the 6 BA stars in our sample range (excluding A6) from [Z] = —0.5 to 
— 1.0 dex. The weighted mean metallicity, is [Z] = —0.87 ± 0.06 dex, where the uncertainty 
is given by the standar d deviation of the sample. This value compares well with the results 
of iBresolin et al.l (120061 ) . based on the a-element content of three early B-type supergiants. 



with weighted mean of [0/H ] = — 0.86 ± 0.0 7 dex. It also agrees with the oxygen abundances 
obtained from H II regions (ILee et al.ll2005l ). This agreement indicates that WLM's young 
stellar population while metal poor exhibits an abundance pattern very similar to the Sun 
(at least with respect to the most relevant species, like O, Fe or Ti). 

It seems that Al4's metallicity is an outlier, being 0.26 dex above the mean, which is 
a bit larger than the individual uncertainties we are claiming in our analyses. As discussed 



above, our low resolution result is consistent with the high resolution analysis by I Venn et al. 



(120031 ). The difference with re spect to the r est of our sample, along with the abnormally high 
oxygen abundance derived by lVenn et al.l (120031 ) could indicate that the star has a peculiar 
history. We note that K08 in their study of NGC300 have found two similar outliers not 
representing the expected metallicity at their galactocentric distance. They expressed that 
the idea of homogeneous metallicity might be naive. Excluding A14 from the weighted mean 
results in only a slight change to [Z] = —0.89 ± 0.07 dex. 



Combining the results of our sample and the early B-type supergiants from lBresolin et al. 



(120061 ) . there is no evidence of a metallicity dependence of the young stellar population with 
the spatial location in the galaxy (see Fig. [TTl) . 
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5.2. Extinction 



Redd ening; values were det ermined u sing Johnson V-band an d Cousins I-band photom- 
etry from iBresolin et al . (see also iPietrzyhski et al.l 120071 ). along with B-V and V-R 
colors from lMassey et al.. (2007) . V- and I-band data for the stars in common in these papers 
are shown in Fig. [71 For our sample of 11 B- and A- type supergiant stars (represented by 
filled dots in the figure) we obtain a difference in the zero points of —0.03 ± 0.02 mag and 
— 0.013±0.030 mag in V and I respectively, with the magnitudes from the first reference being 
fainter in both filters. Note that the star A17 was not considered in the I-band mean since 
it presents a large difference of almost -0.5 mag betwe en both datase t s. Us ing multi-epoch 
photometry (around 100 epochs, spanning over 2 yr), iBresolin et al.l (120061 ) did not detect 
variability, beyond the observational scatter, for any of our targets. It is thus unlikely that 
the discrepancy in A17 I-band magnitudes is related to variability. In the case of the early 
B-type supergiant A9, its B-band magnitude is not consistent with all the other photometric 
measurements (see the corresponding plot in Fig. ITUl) . In both cases, we did not take these 
values into account for the determination of the corresponding extinctions for the stars. 

For four of t he stars in th e samp le (A4, A6, A9, and AlO), we also have IR J- and 
K-band data from lGieren et al.l (120081 ). 



The individual reddening values are presented in Table HI and some examples comparing 
the reddened synthetic SEDs with the photometry are shown in Fig. [TOl Individual reddening 
values range from 0.03 to 0.12 mag, with the extreme case of 0.29 mag for A6. Excluding 
this object, the mean reddening of the BA sample is 0.07 ±0.01 mag, or 0.08 ±0.01 mag 
when considering also the three early B-type supergiant s. This mean value i s higher than the 
characteristic foreground value of E(B-V)=0.037 mag (ISchlegel et al.lll998l ). indicating that 
our stars suffer from inter nal (WLM produced ) reddening. Our mean value is in agreement 
with the recent result by iGieren et al.l (120081 ) based on multi-wavelength observations of 
Cepheids. These authors find a mean characteristic reddening of 0.082 ± 0.02 mag. 

With regard to A6, we note that a very high reddening value is derived for both the 
high temperature and the low temperature solution (discussed in Sect. 13.51) . This means that 
the star suffers three times the mean reddening. A6 could be spatiall y associated with, o r 
close enough to, one of the high column density H I areas identified by lKepley et al.l (120071 ). 
in particular the region identifie d by these auth o rs as the handle of the hook, extending 
north of the CI Ha complex of iHodge &: Millerl (Il995l ). However, two of the three early 
B-type supergiants, A9 and AlO, seem to be also spatially related to high column density 
H I regions, but they do not present such extreme reddening. We cannot discard, however, 
the possibility of a patchy ISM on small scales, and comparable cases of high reddening have 
also been found by K08 in NGC300. 
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5.3. Masses, radii and luminosities 

Once the fundamental parameters Te// and log g are known, it is possible to derive 
the distance dependent quantities by adopting a distance to WLM. In the next section, 
we will derive a distance to the galaxy based on the Flux weighte d Gravity-Luminosity 



Relationship, FGLR, of blue supergiant stars (iKudritzki et al.ll2003l ). For the purposes of 
this section, we adopt that distance, (m — M)q = 24.99 ± 0.10 mag, without any further 
explanation, delaying its derivation to the next section. 

Using this distance, the apparent bolometric magnitudes can be converted to absolute 
bolometric magnitudes, and luminosities. From the luminosities and the effective tempera- 
tures, we compute the radii. With surface gravities and radii, it is then possible to derive 
(spectroscopic) masses. All these quantities are summarized in Table [51 To estimate their 
uncertainties, we propagate the uncertainties in the fundamental parameters. 

From the luminosities, we can also derive stellar masses using theoretical evolutionary 
models (see Fig. [8]). For this purpose, we use the mass-luminosity relatio nships presented by 



K08. In par ticular, we selected K08's fits to SMC metallicity models from lMaeder fc Meynet 



( I2OOII ) and iMeynet fc Maeden (120051 ) that take into account rotation. These evolutionary 
masses are given in the last column of Table O The uncertainties of the spectroscopic 
masses take into account uncertainties in luminosities, gravities and radii. The uncertainties 
in evolutionary masses account only for the uncertainties in the derived luminosities and are 
much smaller. They are, however, affected by the systematic uncertainties of evolutionary 
tracks. Fig. [9^ shows the comparison of evolutionary and spectroscopic masses, including the 
objects in NGC300 studied by K08. Given the typical spectroscopic mass errors of ±0.1 to 
0.2 dex, we conclude that we have agreement between the two types of mass determinations. 
However, there is an indication that the spectroscopic masses are systematically smaller 
than evolutionary masses for high luminosities (above Mboi ~ —8 mag, see Fig. [9)d). We will 
discuss this in section 16.11 

From the direct inspection of the HRD and the masses given in Table we conclude 
that the BA stars of our sample were born with masses in the range of 8 to 20 Mq, with the 
early B supergiants evolved from slightly more massive progenitors between ~25 to 50 Mq. 
This is in agreement with the results found by K08 for NGC 300, who explain the difference 
of the masses as a selection effect. 



6. Flux weighted gravity— luminosity relationship and distance to WLM 

Kudritzki et al.l (l2003l . see also K08) revealed the existence of a tight correlation between 
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the flux weighted gravity, gp = g/T^g, and the luminosity of BA supergiant stars. This 
relationship is supported by the predictions of evolutionary models. Very briefly, the physical 
reason behind this is that during their evolution from the Main Sequence, massive stars with 
masses below ~ 60 Mq will evolve at almost constant luminosity and, due to the short 
timescale and the low mass loss rates, constant mass. In this case, the luminosity of the star 
is correlated with the flux weighted gravity. 

K08 have found a relation between the absolute bolometric magnitude and the flux 
weighted gravity of the form 

Mboi = a(log(?F-1.5) + b (1) 

with their present best values of the coefficients a and b derived from a large sample of su- 
pergiants combining spectroscopic results obtained for 8 different galaxies. These coefficients 
are given in Table [61 We also include the values obtained by K08 for only the supergiants 
in NGC 300, to give an idea of the uncertainties involved in calibrating this relationship. 

The supergiants in WLM show a tight fglr when plotted in apparent bolometric mag- 
nitudes (in Fig. [12]). This can now be used to determine a distance. Following K08, we fit 
an expression of the form 

"^boi = a {loggp - 1-5) + b 

This fit is shown in Fig[T2]by the solid line, and the values for a and b are given in Table [2] 
as well. With only ten stars the slope is somewhat uncertain, thus, we prefer to fix the slope 
to the value obtained by K08 for their large sample. We then recalculate the fit with this 
fixed slope (dashed line; the new value for the coefficient b is also presented in Table [6]). The 
difference in b with respect to K08 then yields a first determination of the distance modulus, 
fi = {m — M)g, for which we obtain /x = 24.99 mag. 

If we fix the slope to K08's fglr for NGC 300 stars only, and compare again, we 
derive /i = 25.06 mag. This can be interpreted as an estimate of the possible systematic 
uncertainties to be about 0.06 mag. 

The statistical uncertainty of the distance modulus is given by 

n(n — 1) ^ af 

^ ' i=i « 

where n is the number of WLM stars, M^^^ is obtained from Eq. [T] evaluated at the 
corresponding observed (7F,i, ^boi,i is the apparent bolometric magnitude vnboi^i corrected by 
the derived distance modulus /x, and di = ai/ Ylij ^j^ with a} = (^m^^i'^^'^^ioggp individual 
uncertainties in mf,oi, accounting for the uncertainty in the BC, in the extinction and in the 
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observed photometric errors, as well as the uncertainties in the flux weighted gravity. We 
obtain a statistical uncertainty associated with the distance modulus of 0.10 mag. 

An alternative way to determine the distance modulus is to minimize the residuals in 
magnitudes once the stars are shifted to a particular distance, i.e. we determine the value 
of /i that minimizes 




For this second method, the uncertainty in the distance determination is given by the square 
root of the minimized residuals. 

The distances derived by both methods are in excellent agreement. This is a consequence 
of the fact that the slope of the relationship defined by the WLM stars is very close to the 
slopes found by K08 (see Tabl^. Fig. [T3l presents the final distance corrected fglr for 
the WLM stars (filled circles), together with all the objects used by K08 to calibrate the 
relationship. This figure also includes K08's fglr fit used to calculate the distance to WLM. 
As can be seen in this figure, the agreement is very good. 

Finally, we would like to point out that the solution used for A6 (cool or hot model) 
has little effect on the derived distance, since the values of the flux weighted gravity are very 
similar in both cases, and in any case, a change in the flux weighted gravity is correspondingly 
accompanied by changes in the bolometric correction and extinction (relatively minor in this 
case), in such a way that the star moves along the relationship. The physical reason behind 
this behavior of gp has been extensively discussed by K08, to which the reader is referred 
for further insights. 



6.1. WLM FGLR: empirical versus theoretical relationships and metallicity 

dependence 



In Fig. [H] we show the observed fglr for WLM compared with the prediction of stel- 
lar evolution. We have chos en the stellar evolu t ion f glrs for solar and SMC metal licity 
using evolutionary track s by iMeynet fc Maederl (120031 ) , and iMaeder fc Meynetl (120011 ) and 
Meynet fc Maederl ((20051), which include the effects of stellar rotation with initial rotational 
velocities of 300 km s~^ (for the related mass-luminosity relationships and the parameteriza- 
tion of the evolutionary FGLRs see K08). 



We note that for low luminosity and high flux-weighted gravity [log gp > 1.6, Mboi < 
—8 mag) there is good agreement. However, towards lower gravities and higher luminosities 
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the stellar evolution FGLRs show a strong curvature and start to differ from the result of 
the spectral analysis. The effect is stronger for SMC metallicity ([Z]=-0.7), which is close 
to the average metallicity we have determined for the young stellar population in WLM. An 
indication of a similar discrepancy has already been noted by K08 (see their Fig. 27 and the 
corresponding discussion). 

This discrepancy seems to be equivalent to a discrepancy between stellar evolutionary 
and spectroscopic mass at high luminosities (see Fig. [H]). Note that (see K08, Eq.(29)) 

(7F = oc MR-'T-^' oc ML-' 

and, in consequence, smaller spectroscopic masses result in a fglr shifted towards smaller 
flux weighted gravities. Thus, one reason for the discrepancy could be a systematic un- 
derestimate of stellar gravities by the spectral analysis at the high luminosity end of the 
FGLR. We realize that the most discrepant objects in Fig. [9] and Fig. [H] are early B-type 
supergiants, which where analyzed with a different stellar atmosphere code. On the other 
hand, K08 had a much larger sample of supergiants available and did not find a systematic 
difference between high luminosity early B-types and BA-types. 

Another possibility might be the effects of rotation on the evolutionary tracks. K08 
noticed that tracks without rotation result in FGLRs shifted towards higher gravities because 
the mass-to-light ratio is larger. Thus, tracks with even higher initial rotation than consid- 
ered here may result in FGLRs in better agreement with the observations. However, such 
increase in rotational velocity will not be supported by recent studies of LMC and SMC 
Main Sequence B stars, progenitors of BA supergiants, sho wing projected rotational veloci- 



ties below 200 km s ^ (see for example iHunter et al.ll2008l ) 



Finally, it is also important to discuss a possible metallicity dependence of the fglr. 
Except for the very high luminosity end, the relative difference between the evolutionary 
FGLRs for [Z]=0.0 and [Z]=-0.7 is very small. This difference is mostly caused by the effects 
of mass-loss, which are stronger at higher metallicity and higher luminosity. The fact that 
we obtain a distance consistent with the TRGB and Cepheid studies of WLM using the 
FGLR calibration by K08, which is based on objects mostly with LMC metallicity, supports 
the conclusion that metallicity effects are not important. Future work on supergiants in 
other metal poor galaxies will very likely help to clarify the situation. 
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7. Conclusions 

The primary goal of this work has been to determine the distance to WLM using a new 
spectroscopic method, the fglr. The distance modulus obtained, 24.99=fc0.10 mag, com- 



pares well with recent determinations, based on purely photometric methods. iRizzi et al. 



( 120071 ) found 24.93±0.04 mag from the TRGB (see also iPietrzyhski et al.l 120071. an d refer- 
ences in Sect. [1]) and the multi-wavelength study of Cepheids by lGieren et al.l (120081 ) yielded 
24.924±0.04±0.04 mag (statistical and systematic errors). With only 10 stars available for 
our study, the statistical uncertainty is larger than the ones claimed by the photometric 
studies. However, this could certainly be improved with a larger sample of objects. Our 
study also confirms that the young population of WLM suffers from significant extinction, 
higher than the foreground value. This is another example where the accurate measurement 
of intrinsic extinction turns out to be important for the determination of distances. It is the 
advantage of the spectroscopic method presented here that it yields reddening and extinction 
for the individual targets. 

We have also determined stellar metallicities and found an average value of [Z]=- 
0.87 ±0.06 dex, in good agreement with previous H II region studies. The fact that the 
TRGB and fglr distances agree at this low metallicity indicates that metallicity effects are 
very likely small for the fglr method. 

Stellar evolution calculations can be used to determine a theoretical fglr, which can 
then be compared with observations. We find good agreement at low luminosities indepen- 
dent of metallicity assumptions for the theoretical relationship, but disagreement at higher 
luminosities with the theoretical fglr, which are more important for the low metallicity. 
Assuming higher initial rotational velocities, which would enhance mass-loss and rotational 
mixing, might be a way to solve this discrepancy. These higher rotational rates are not 
supported by recent studies of the progenitors of these objects in the Magellanic Clouds. 

Finally , our result s for A 14 (aka WLM-15) are in good agreement with the results 



obtained by lVenn et al.l (120031 ). confirming to some extent the particular nature of this object. 



Unfortunately, we could not perform an independent analysis of its oxygen abundance. 

In summary, this first application of the FGLR method to determine a distance to a 
galaxy has proven to be successful. The fglr technique seems to be a robust and reliable 
way to provide independent and accurate information about extragalactic distances. 
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A. Synthetic photometry 



In this appendix we present specific details relative to our synthetic photometry, that 
we consider could be useful to others. 



We follow the ideas presented by iBesselll (l2005l . in particular see its Section 1.6) in 
order to compute the photometric magnitudes from our model atmosphere models. Given a 
spectral energy distribution fx, the magnitude in a given filter with transmission curve 
Ra (A) is calculated accounting for the fact that modern detectors count the number of 

photons, not energy, therefore tua oc fx {XRa (A)) dX. 



In order to determine the zeropoints of the differen t band passes, we use the recent spec- 
tral energy distribution of Vega presented by iBohlinI (120071 ). along with the filter curves 
obtained from different references. Table [7] contains the zeropoints of our synthetic photom- 
etry, as well as the references for the filter transmission curves. 
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Table 1 . WL M stars analyzed in this work. Primary identifica tion following 



Bresolin et al.l (120061 ) . Alternative identification and photometry from lMassey et al.l (120071 ) 
is also included. For c omple teness, information about the three early B supergiants studied 
by lBresolin et al.l (120061 ) is presented in this table. All photometric quantities are in 

magnitude units. 



ID 


V 


V-I 


Spectral Type 


S/N 


Alt. ID 


V 


B-V 


V-R 


R-I 




(Bresolin ct al. 2006) 




(Massey et al. 2007) 






A6 


19.82 


0.38 


A7Ib 


49 


0156.16-152624.5 


19.789 


0.217 


0.179 


0.203 


A4 


20.22 


0.07 


A2II 


44 


0201.57-152527.0 


20.185 


0.016 


0.038 


0.056 


A14 


18.43 


0.23 


A2II 


96 


0159.56-152926.1 


18.374 


0.087 


0.073 


0.120 


A16 


18.44 


0.16 


A2 la 


96 


0157.89-153013.1 


18.383 


0.204 


0.058 


0.055 


k2^ 


20.16 


0.09 


AO II 


44 


0159.04-152442.8 


20.142 


-0.014 


0.052 


0.027 


A12 


17.98 


0.06 


B9Ia 


119 


0153.22-152839.5 


17.966 


0.005 


0.031 


0.014 


A5 


19.41 


-0.04 


B8Iab 


64 


0203.31-152552.6 


19.412 


-0.117 


-0.021 


-0.082 


A17 


19.34 


0.00 


B5Ib 


67 


0200.81-153024.8 


19.313 


-0.109 


-0.026 


0.451 


Agi^ 


18.44 


-0.06 


B1.5Ia 


101 


0157.20-152718.0 


18.392 


0.201 


-0.028 


-0.059 


AlO'" 


19.34 


-0.15 


BO lab 


68 


0154.06-152745.4 


19.317 


-0.171 


-0.050 


-0.093 


All^ 


18.40 


-0.18 


09.7Ia 


106 


0159.95-152819.0 


18.378 


-0.109 


-0.069 


-0.120 



^This stars is outside the limits of the grid of models, therefore it cannot be analyzed. 
''Early B-type supergiant studied by lBresolin et a~ I ll2006h . 



Table 2. Mean J and K photometry from multi-epoch observations by lGieren et al.l (120081 ) 



ID 


J 


a (J) 


K 


a{K) 




(mag) 


(mag) 


(mag) 


(mag) 


A6 


19.055 


0.022 


18.896 


0.129 


A4 


19.889 


0.021 


19.843 


0.067 


A2 


19.916 


0.010 


19.748 


0.158 


A9 


18.578 


0.015 


18.709 


0.111 


AlO 


19.624 


0.009 


19.673 


0.082 
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Table 3. Fitness of the solutions. See text for definitions of it! and Q. Stars are grouped in 
two sets, according to the method used to determine the temperature. For the first group, 
the Balmer jump is used, and for the second we use the He I hnes for A5 and A17, and the 

relative strength of Ti II hnes for A4 (see text) . 



ID 


R (log g) 


R{[Z]) 


R {Teff) 


R 


S/N 


Q {log g) 






A6^ 


0.1273 


0.1397 


0.0799 


0.3469 


49 


-0.3480 


4.2911 


-0.9223 




0.1190 


0.1435 


0.0213 


0.2838 


49 


1.0856 


4.3364 


-0.0874 


A14 


0.1085 


0.0430 


0.0210 


0.1726 


96 


-0.6183 


0.2116 


-0.0461 


A16 


0.0651 


0.0535 


0.0037 


0.1223 


96 


-0.3119 


0.3723 


0.0067 


A12 


0.0397 


0.0388 


0.0196 


0.0982 


119 


-1.0653 


0.0388 


0.0675 


A4 


0.1404 


0.1047 


0.0836 


0.3288 


44 


0.2761 


-1.7217 


-0.0375 


A5 


0.0641 


0.0768 


0.0689 


0.2099 


64 


-0.7499 


0.8870 


0.0169 


A17 


0.0842 


0.0479 


0.0305 


0.1682 


67 


-0.5869 


0.2315 


0.0004 



"Cool solution, does not reproduce the Balmer jump, see text 
''Hot solution that reproduces the Balmer jump 
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Table 4. WLM - Stellar parameters 



ID 




log g 


log g-e 


[Z] 


E(B-V) 




comments 




(K) 


(cgs) 


(cgs) 




(mag) 


(km s"-"-) 




A6 


8750±240 


1.95±0.17 


2.19±0.12 


-0.15±0.33 


0.295±0.010 


4 


a 


A4 


8550±350 


2.00±0.20 


2.27±0.13 


-0.77±0.26 


0.060±0.030 


4 


b 


A14 


8270±144 


1.60±0.12 


1.93±0.09 


-0.50±0.19 


0.124±0.007 


4 




A16 


10650±120 


1.78±0.05 


1.67±0.03 


-0.70±0.23 


0.112±0.023 


6 




A12 


12100±136 


1.79±0.03 


1.46±0.01 


-0.78±0.21 


0.061±0.004 


8 




A5 


12220±430 


2.20±0.07 


1.87±0.02 


-0.70±0.14 


0.030±0.010 


4 


c 


A17 


13500±450 


2.35±0.07 


1.86±0.03 


-1.00±0.15 


0.061±0.004 


5 


c 


A9 


20000±1000 


2.45±0.10 


1.62±0.03 


-1.00±0.20 


0.130±0.020 


12 


d 


AlO 


25000±1000 


2.90±0.10 


1.62±0.03 


-0.80±0.20 


0.160±0.020 


15 


d 


All 


29000±1000 


3.00±0.10 


1.62±0.03 


-0.80±0.20 


0.070±0.019 


15 


d 



"Metallicity uncertain, see text 

''no Balmer jump measured; T^f f determined from spectrum (Ti II lines, see text) 
■^no Balmer jump measured; Tgff from Hel lines, see text 

'^Results from [ Bresolin et al] ll2006h . with E(B-V) updated with the IR photometry as well as 
iMassev et al. I bom data. 

Note. — Note that the microturbulence is not derived in this work. See text for an explanation. 
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Table 5. Bolometric corrections and distance dependent magnitudes: radii, luminosities 

and masses 



name 


mv 


BC 




log L/Lq 


R 








(mag) 


(mag) 


(mag) 


(cgs) b 




(R©)'' 


(M0)'> 


A6 


19.82±0.03 


-0.02±0.02 


18.85±0.04 


4.35±0.05 


65.5±5.4 


14.3±6.1 


10.2±0.4 


A4 


20.22±0.05 


-0.05±0.05 


19.98±0.12 


3.90±0.09 


40.7±5.3 


6.1±3.2 


7.7±0.4 


A14 


18.43±0.02 


O.OOitO.Ol 


18.05±0.03 


4.67±0.05 


106.3±7.2 


16.4±5.1 


12.8±0.5 


A16 


18.44±0.02 


-0.43±0.01 


17.65±0.07 


4.83±0.07 


76.8±6.4 


13.0±2.6 


14.4±0.8 


A12 


17.98±0.01 


-0.70±0.01 


17.07±0.03 


5.06±0.05 


77.8±5.1 


13.6±2.0 


17.4±0.8 


A5 


19.41±0.02 


-0.81±0.03 


18.59±0.08 


4.46±0.07 


38.1±4.1 


8.4±2.3 


ll.OitO.5 


A17 


19.34±0.02 


-0.85±0.02 


18.28±0.07 


4.58±0.07 


38.9±4.1 


11.3±3.0 


11.9±0.6 


A9 


18.44±0.02 


-1.91±0.12 


16.13±0.18 


5.44±0.11 


43.9±7.1 


19.8±7.9 


26.8±2.7 


AlO 


19.34±0.02 


-2.46±0.10 


16.39±0.21 


5.34±0.12 


25.0±4.0 


18.1±9.2 


24.5±2.7 


All 


18.40±0.02 


-2.81±0.09 


15.38±0.17 


5.74±0.11 


29.5±4.0 


31.8±11.6 


35.9±4.0 



"Apparent bolometric magnitude: mdoi = - ^v) + BC 

''Distance dependent magnitudes evaluated for (m — M)q = 24.99 ± 0.10 mag 
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Table 6. FGLR coefficients, {Mboijmboi} = o, {log g-p — 1.5) + b. 



Relationship 


a 


b 


Comments 


K08-all 


3.41±0.16 


-8.02±0.04 


Mboi, calibrating fglr using stars in 8 galaxies 


K08-NGC 300 


3.52±0.25 


-8.11±0.07 


Mboi, calibrating FGLR using only NGC300 stars 


WLM 


3.48±0.36 


16.95±0.12 


™bol) WLM stars analyzed in this work 


WLM-all 


3.41 


16.97±0.09 


™bol) WLM stars with FGLR slope from K08-all 


WLM-NGC 300 


3.52 


16.95±0.08 


mbol, WLM stars with fglr slope from K08-NGC300 



Note. — 



Adopted slopes are shown in italics 
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Table 7. Synthetic photometry zeropoints 



Bandpass 


Zeropoint 


Ao 


Reference 




(mjy) 


(A) 




V 


3647.62 


5450 


Bessell (1990) 


I 


2432.91 


7980 


Bessell (1990) 


B (KPNO) 


4070.71 


4381 


NOAO filter K1002 


V (KPNO) 


3726.61 


5387 


NOAO filter K1003 


R (KPNO) 


3100.14 


6513 


NOAO filter K1004 ^ 


I (KPNO) 


2488.43 


8245 


NOAO filter K1005 ^ 


J (2MASS) 


1594.00 


12350 


2MASS web page ^ 


Ks (2MASS) 


667.00 


21590 


2MASS web page ^ 



''http: / / www.lsstmail.org/kpno/mosaic / filters / filtcrs.html 
''http: / / spider.ipac.caltech.edu / staff/waw / 2mass / opt_cal / 
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4850 



4860 



4870 




4330 4340 4350 




4090 4100 4110 
Wavelength (A) 



o 



2.80 




+ 






0.08 


2.90 




+ 

+ 

+ 






0.06 


3.00 
3.10 


4 

+ 

+ 

+ 


+ 




;(S-0)I 


0.04 


3.20 








0.02 


3.30 


+ 


1.46 dex 






0.00 
1 


1 


1 1.2 1.3 


1.4 1.5 1.6 1 


7 






1.46 dex 
1 1.2 1.3 1.4 1.5 1.6 1 



2.70 






+ 

+ 






0.08 


2.80 




+ 

+ 








0.06 


2.90 




• 

+ 






m 






+ 








o_ 


0.04 


3.00 


+ 














+ 

+ 












3.10 










0.02 




+ 












3.20 




1.45 


dex 






0.00 


1 

2.60 


1 1.2 1.3 


1.4 1.5 


.6 1 


7 




1 












0.06 


2.70 






+ 

+ 






0.05 


2.80 




+ 

+ 








0.04 


2.90 




+ 


















0.03 




+ 












3.00 


+ 










0.02 


3.10 


+ 

+ 










0.01 




+ 












3.20 




1.41 


dex 






0.00 


1 


1 1.2 1.3 


1.4 1.5 


.6 1 


7 




1 




log 


g (dex) 












1.46 dex 



1.1 1.2 1.3 1.4 1.5 1.6 1.7 




1.41 dex 



1.2 1.3 1.4 1.5 1.6 1.7 
log g (dex) 



Fig. 1. — Determination of surface gravity. For illustration purposes, only three lines are 
presented. From top to bottom: H7 and US. For each line, the first panel shows the 
base vector, the second one displays the projection of the models (crosses) and the observed 
spectrum (filled dot), and the third panel presents a solution by applying a minimum distance 
method (minimization of the quadratic differences, (Oa — S{) ) for comparison. The 
derived log g for each line is also included in each individual plot. 
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Fig. 2. — Metallicity determination. For illustration purposes, only four spectral windows 
are shown. For details about each individual panel, see previous figure. 
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Fig. 3. — Determination of T^ff-. (a) definition of the Balmer jump region; diamonds identi- 
fied the wavelength points use to create the PCA base vector, (b) Tgff determination based 
on the minimum difference between the observed and modeled Balmer jump measurements, 
with Db = < logPf" > - < logFf > following K08; the cross marks the derived temper- 
ature, (c) PCA base vector for the Balmer jump, and (d) projection of the models (hollow 
diamonds) and the observed Balmer area (cross). 
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Fig. 4. — Teff-log g diagram for star A14 of our sample. The solution provided by the 
algorithm is defined by the small box, while the solution obtained by the conventional method 
would be located at the intersection of the two T^jf-log g fit curves for the Balmer jump 
(x) and the Balmer lines (+). The dashed lines define the final uncertainties in the derived 
parameters. 
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Fig. 5. — Relevance of the S/N on metallicity determination. First four plots correspond to 
a [Z] = 0.00 dex case; the second set of plots present the case for [Z] = —0.85 dex. In each 
individual plot, the x-axis represents the difference between input and derived metallicities, 
and the y-axis gives the relative frequency, for a total number of 100 independent trials. A 
Gaussian fit to the resultant distribution is shown, and its sigma is given in each plot. 
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Fig. 6.— UVES/VLT (R=32000) spectrum of WLM-A14 compared with our solution (red) 
obtained from the analysis of low-resolution F0RS2/VLT (R~1000) spectrum. 
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Fig. 7. — Comparison of photo metri c measuremen t s. Fir st row: V- and I-band photometric 
data from lBresohn et al.l (120061 ) and lMassey et aLl (120071 . M2007 subindex), for a total num- 
ber of 35 stars in cor nmon (diamonds'). Th e stars analyzed in this work, and the three early 
B-type supergiants of iBresolin et al.l (12006! ) are represented by filled circles. The 1:1 relation- 
ship is represented by the dashed line. Second row: considering only the 11 stars analyzed 
in our work, we find a mean difference in the zero point calibration of —0.03 ± 0.02 mag and 
—0.013 ± 0.030 mag in V and I respectively. The star A17, not considered in the I-band 
mean, presents a difference with respect to our reference value of almost -0.5 mag. 



-40- 



G 



6.0 
5.5 

5.0 



6OM0 



40M 



° 4.5 



4.0 
5.5 



4.1 




O.2Z0 (Maeder & Meynet 200i; 



4.6 4.4 4.2 4.0 5.8 5.6 
log Teff 



Fig. 8. — WLM Hert zsprung-Russel diagra m. Circles locate the three early B-type su- 
pergiants analyzed by iBresolin et al.l (120061 ). once corrected for the difference in distance 
modulus. The filled squares display the sam ple of BA supergiants ari alyzed in this work. 
Evolutionary tracks for rotating models from iMaeder fc Meynetl (120011 ) are also shown, la- 
beled with their initial mass. 
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Fig. 9. — WLM stellar masses, (a) Comparison of evolutionary a nd spectroscopic mass es for 
our sample of WLM stars (filled circles) and the B supergiants of iBresolin et al.l (120061, filled 
squares). The sam ple of NGC300 BA su pergiants from K08 (circles) as well as the early 
B supergiants from lUrbaneja et al.l (l2005l . empty squares) are also shown, with the dotted 
line defining the 1:1 relation, (b) The ratio of spectroscopic to evolutionary mass versus the 
stellar luminosity, for the WLM sample. Symbols as in (a). 
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Fig. 10. — Synthetic reddened SEDs and photometric measurements for different stars 
in our WLM c ombin ed sample. Circles correspond to V and I photo r netric values from 
Bresolin et al.l (|2006[) a s well as J and Ks photometry from iGieren et al.l (120081 ) , and crosses 



stand for 



Massev et al.l (120071 ) B,V,R and I-band data. 
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Fig. 11. — Spatial distribution of the objects. The image is a false color composition us- 
ing [S II] (red channel), Ha (green channel) and [O I II] (blue channel) narrow band im- 
ages from the Local Group Survey (jMassey et al.l 120071 ). Star's color coding: white circles 
E(B-V)<0.05, cyan 0.05<E(B-V)<0.10, yellow 0.10<E(B-V)<0.15, and white boxes E(B- 
V)>0.15 mag. North is up and East is to the left of the image. 
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log (cgs) 

Fig. 12. — WLM FGLR based on apparent bolometric magnitudes, mboi = iti^ — + BC. 
The solid line is the best linear fit to the data, while the dashed line shows a linear fit when 
the slope is fixed to K08's FGLR based on stars in 8 different galaxies (see text). The star 
A6 is identified in this figure by the square. 



-45 - 




Fig. 13. — WLM FGLR. This figure presents the bolometric magnitudes of our WLM stars, 
once corrected for the distance modulus derived through the FGLR. The sample of K08 is 

represented by open symbols (circles: NGC 300 stars; diamonds: other galaxies. See K08 for 
details). The thick line corresponds to the FGLR calibration used to determine the distance 
to WLM. 
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Fig. 14. — WLM FGLR. Comparison of empirical and theoretical relationships. The 
theore tical fglrs are obtaine d from evolutionary mod els with rotation for SMC metal- 
l icitv f|Maeder fc Mevnet 2001 : Meynet Sz Maeder 2005 . dashed line) and solar metallicity 
( Meynet fc Maeder 20031 . dot-dash line). See text for a detailed discussion. 



